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Abstract 



Heavy quarkonium production is expected to be sensitive to the formation of a 
quark gluon plasma (QGP). It was (and still is with ongoing data analyses) exten- 
. sively studied at the CERN SPS, at collision energy ^/snn of the order of 20 GeV. 

An anomalous suppression was clearly observed. The PHENIX experiment at RHIC 
has presented preliminary results that exhibit a similar amount of J/tp suppression, 
at ten times higher collision energy. I review the results obtained at both facilities. 
While interpreting and comparing them, the importance of understanding normal 
nuclear effects is emphasized. A new method to derive a reference for Au+Au colli- 
sions from the centrality dependence of d+Au measurements at RHIC is presented. 
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^ . 1 Nuclear effects from p+A or d+A collisions 

.W. 

When trying to interpret all the quarkonia yields observed in p+A collisions 
at various energies and kinematical domains, we find ourselves facing a real 
puzzle [1]. Various effects are invoked, including normal nuclear absorption 
of quarkonia or pre-resonant cc pairs, parton shadowing and corresponding 
anti-shadowing, contribution from the intrinsic charm existing in the nucleon 
wave function, or energy loss and related transverse momentum broadening, 
all of them being further complicated by feed-down from higher mass states. 
Explaining all the available p+A data is beyond the scope of this article. In 
the following, I restrict my interest for p+A data and nuclear effects to the 
energies and kinematical domains where A+A collisions are also measured, in 
order to set up references for QGP studies. 



Email address: raphael@in2p3.fr (Raphael Granier de Cassagnac) 



Preprint submitted to Elsevier Science 



8 February 2008 



i 80 

>r to- 

8 60 
° 50 

I 40 



SQ " 



30 



20 



10 - 

9 



Experimental rescale to 200 GeV 
Theoretical rescale to 158 GeV 

NA51 pp, pd 450 GeV 
□ NA50 LI p-Be, Al, Cu, Ag, W 450 GeV 
NA50 HI p-Be, Al, Cu, Ag, W 450 GeV 
A NA50 p-Be, Al, Cu, Ag, W, Pb 400 GeV 




NA38 S-U 20(1 GeV 



T NA50 Pb-Pb 2000 158 GeV 
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Fig. 1. Normal nuclear matter effects. Left) At SPS, the J/^/Drell-Yan cross-section 
ratio versus the average nuclear length L, for several collision systems [2]. The line 
stands for a normal nuclear absorption of ~ 4.2 mb. Right) At RHIC, J/ijj nuclear 
modification factor RdA as a function of the number of nucleon-nucleon collisions 
N co u for backward (top, y = —1.7) mid (middle, y = 0) and forward (bottom, 
y = 1.8) rapidities [6]. Theoretical curves from Vogt [7], assuming EKS shadowing 
and 0, 1, 2 or 3 mb normal absorption cross-sections (from top solid to bottom 
dot-dashed). 

1.1 Normal nuclear absorption at SPS 



Figure 1 (left) summarizes the published J/ip production yield (normalized to 
Dreil-Yan) observed at SPS energies, from p+p to Pb+Pb collisions [2], exclud- 
ing In+In data from the NA60 experiment [3]. It is plotted here as a function 
of the average length L of nuclear matter traversed by the cc state. Looking 
only at p+A collisions, we observe a clear exponential suppression. This be- 
havior is expected if the only nuclear effect is an absorption by the nucleons in 
the incoming nuclei. The line of figure 1 (left) stands for such a normal nuclear 
absorption, with a fitted J/^-nucleon cross-section of a 3 J^ = 4.18 ± 0.35 mb. 
Although this simple picture does a splendid job in reproducing the data, all 
the way from p+p up to peripheral Pb+Pb collisions (including the whole S+U 
range), one can oppose theoretical arguments against it. First, the p+A data 
are rescaled to a lower collision energy ^s^n (from 27 or 29 GeV to 17 and 
19 GeV). The NA60 collaboration is currently analyzing p+A data at 17 GeV. 
Second, according to various theoretical predictions [4], shadowing (or rather 
anti-shadowing) could play a role at the SPS experiments regime, namely for 
momentum fraction x of the order of 10 _1 . The SPS experiments can hardly 
address this question because of their limited rapidity (or x) range. Still, the 
rapidity asymmetry observed in J / ip yields in p+A collisions might be a hint of 
such an effect [5] . Third, the p+A absorption might not be straightly applica- 
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ble to A+A, for instance if there are complications due to changing feed-down 
ratios. Indeed, the ip 1 is known not to suffer the same nuclear effects as the 
J/ip: its absorption cross-section extracted from p+A is a^ bs = 7.7 ± 0.9 mb 
and its anomalous suppression already sets in S+U collisions [5]. However, 
one should not forget that the simple picture of nuclear absorption perfectly 
describes a large amount of data (p+p, various p+A, S+U and peripheral 
Pb+Pb). The departure from the absorption curve for more central Pb+Pb 
collisions is clearly a new phenomenon. 

1.2 Shadowing and absorption at RHIC 

At ^snn = 200 GeV, the only available data on quarkonia from p+A like col- 
lisions are J/ip seen by the PHENIX experiment in d+Au [6], probing a wide 
rapidity range, from —2.2 to 2.4, corresponding to momentum fractions x of 
gluons in the gold nuclei ranging from ~ 10~ 3 to ~ 10 _1 . The minimum bias 
points of figure 1 (right) show that the measured nuclear modification factors 
RdA depend on rapidity. This is interpreted as due to shadowing and/or anti- 
shadowing. While the strength of gluon shadowing is not heavily constrained 
by theory (models predictions [4] differ by a factor of three), PHENIX data fa- 
vor moderate shadowing schemes such as the Eskola-Kolhinen-Salgado (EKS). 
In addition to this, a moderate normal nuclear absorption is allowed, not larger 
than 3 mb. The addition of these two ingredients is performed by Vogt in [7] 
and can describe both the rapidity and centrality dependencies, as shown on 
figure 1 (right). 

1.3 From d+Au to Au+Au at RHIC 

To interpret J/ip production in A+A collisions at RHIC, we need a model 
capable of reproducing the d+Au nuclear modification factors RdA and ex- 
trapolating them to A+A. Such a model, including inhomogeneous shadowing 
and nuclear absorption is given by [8] and shown on figure 1 (right). An- 
other attempt to derive a reference from d+Au can be found in [9] where the 
authors derive dissociation cross-sections from the centrality dependence of 
RdA- They assume that nuclear effects are proportional to exp(— poadissL), 
Po being the normal nuclear density and L the average length of nuclear 
matter seen by the J/ip. However, there is no fundamental reason for this 
function to reflect the centrality dependence of shadowing. I propose here 
an alternate method, with a concern to be as much data-driven as possible. 
First, I perform phenomenological fits of the modification factors RdA{y,b) 
as a function of the impact parameter b (given by a Glauber model), for 
the three rapidities y of the PHENIX measurements. Given the experimen- 
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tal uncertainties, linear fits are sufficient to describe the data. I then run a 
A+A Glauber model. For each A+A collision occurring at a given impact 
parameter b A A, the positions of the N cM elementary nucleon+nucleon colli- 
sions are randomly distributed (following the nuclear densities) providing the 
locations b\ and b\ of each collision i, relative to the center of nucleus 1 and 
nucleus 2. For the considered A+A collision, the predicted nuclear modifica- 
tion factor is given by the following summation over the elementary collisions: 
RAA{\y\,b AA ) = Y,^=f(RdA(-y,b\)xR dA (+y,b t 2 ))/N coU . This formula assumes 
that a J ftp produced in a A+A collision at a given rapidity y suffers the prod- 
uct of the nuclear effects that were observed in d+A at this rapidity, by the 
ones of the opposite rapidity 1 (equivalent to a A+d collision). This assump- 
tion is correct for the only two effects considered so far to explain RHIC data, 
namely shadowing and nuclear absorption. Quarkonia production is propor- 
tional to the parton distribution functions (pdf) in each nucleus, while the 
average length is the sum of the length in each nucleus, so that the produc- 
tion is finally proportional to pdfi x pdf 2 x exp(— po{L\ + L2)). This method 
has two advantages 2 . First, the statistical and systematical uncertainties of 
the d+A measurements can be directly propagated to the A+A prediction. 
Second, since it is based on a Glauber calculation, it is easy to predict modifi- 
cation factors for experimental centrality classes. This is done on figure 2 where 
predictions are given for the Au+Au PHENIX centrality classes. No system- 
atic uncertainties from the method itself (Glauber parameters) are calculated. 
The amount of predicted suppression is compatible with Vogt's predictions. 
In the forward rapidity case (left) its uncertainty is smaller than the allowed 
variation between the 1 and 3 mb absorption cross-sections and seem to favor 
intermediate ones. In the mid-rapidity case (right) the uncertainty is larger. 
At both rapidities, the J/ ip suppression seen by PHENIX [10] in the most cen- 
tral data is larger than the one predicted by Vogt or by the model described 
above, pointing out that there is an anomalous suppression at RHIC energy. 



2 Anomalous suppressions 

In both SPS (figure 1 left) and RHIC (figure 2) data, the most central A+A 
J I ip measurements depart from the nuclear effect predictions, suggesting that 
other mechanisms are involved. Such an anomalous suppression was early 
predicted by Matsui and Satz as a signature of the QGP [11]. 



1 A similar assumption is made in [9] while summing dissociation cross-sections. 

2 Despite these advantages it will be difficult to apply it at LHC where p+A and 
A+A collisions are planed to be measured at difference collision energies and J ftp 
in ALICE will only be measured at positive or negative rapidity for p+A. 
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Fig. 2. Nuclear modification factor Raa as a function of centrality (given here by 
the number of participants N part ) for Au+Au at y = 1.7 (left) and y = (right). 
Squares are preliminary data from PHENIX [10]. Theoretical curves are nuclear 
effects predictions from Vogt [8] , solid and dashed lines being for 1 and 3 mb normal 
nuclear absorption cross sections, respectively. The circles within the shaded bands 
show the prediction from the model presented here. 

2.1 From SPS to RHIC 



Various models account for the J ftp anomalous suppression seen in Pb+Pb 
collisions at SPS. They were used to predict the expected J/ip yield at RHIC 
energy. Three of these predictions are compared to PHENIX Au+Au prelim- 
inary data [10] on figure 3 (left). In [12] (solid line), J/ift's are absorbed by 
comoving particles (of undetermined partonic/hadronic nature). In [13], the 
authors describe the dynamical interplay between suppression and regenera- 
tion of J/^'s in a QGP. The suppression mechanism is dominant for NA50 
energies and is the only one presented here as a dashed line (see figure 4 for the 
full prediction). In [14] (dot-dashed line), a QGP statistical charm coalescence 
model is used. All three models fail to reproduce PHENIX data, overesti- 
mating the measured suppression. Other models such as percolation [15] also 
over-predicts the suppression, suggesting that new mechanisms take place at 
RHIC energy. It is interesting to note that the same models also failed to 
reproduce the In+In data shown by the NA60 experiment [3]. 



2.2 Alternate explanations for RHIC suppression 



Three classes of models exist that can accommodate the amount of anomalous 
suppression seen in the most central collisions. The accuracy of present data 
and the nuclear effects uncertainty, do not allow to favor one over the other. 
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Fig. 3. RHIC Nuclear modification factor Raa as already presented on figure 2 
versus predictions derived from models reproducing SPS data. Theoretical curves 
are predictions from models that describe the SPS anomalous suppression. Left) 
Most of the models over-predict the suppression [12,13,14]. Right) One model [16] 
reproduces the most central suppression (solid/dashed are without /with plasma, 
top/bottom assume 1/3 mb normal absorption cross-section). 

• Detailed transport: One paper [16], simulating J/ip transport in a hydrody- 
namical model, predicts an amount of suppression that matches the most 
central data. It is shown on figure 3 (right) where the authors have added 
nuclear matter effects (nuclear absorption only, 1 or 3 mb) with respect to 
the published paper. The suppression they obtain is not large probably be- 
cause of their description of the boundary between the QGP and the nuclear 
phase. 

• Sequential melting: An important fraction (30 to 40%) of J/ip's comes from 
decays of excited states (if)', Xc) as if is shown by the HERA-B experi- 
ment [17]. They are taken into account in most of the approaches. Recent 
lattice computations indicate that J/^'s could melt at a much higher tem- 
perature than the one that was originally thought. One possible hypothesis, 
defended in [9], is that, both at SPS and RHIC, only the excited states melt, 
leaving all the initially produced J/ip's untouched. 

• Recombination: At RHIC energies, multiple cc pairs are produced, 10 to 20 
in central collisions [18]. Quark mobility in a deconfmed medium could al- 
low uncorrelated charm quarks to recombine when the QGP fireball freezes, 
raising the quarkonium yield with centrality. A balance between suppres- 
sion and enhancement could lead to the intermediate suppression observed 
at RHIC. Figure 4 (left) shows a collection of predictions from various re- 
combination or coalescence models [13,19,20,21]. Unfortunately, the lack 
of knowledge concerning yields and distributions of the initially produced 
charm quarks, as well as of the recombination mechanism, make these pre- 
dictions hardly predictive. A way to search for recombination is to look at 
its impact on the distributions of kinematical variables, such as transverse 
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Thews et al., nucl-th/0505055 
Grandchamp et al. hep-ph/0306077 (total) 
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Fig. 4. Recombination scenarios at RHIC. Left) Nuclear modification factor Raa as 
already presented on figure 2. Theoretical curves come from various recombination 
models [13,19,20,21]. Right) Mean squared transverse momentum versus N co u, from 
p+p (circles) to d+Au (triangles) and Au+Au (squares) for forward/mid rapidities 
(open/full symbols). Upper/lower shaded bands from [21] stand for direct/fully 
recombined J/tp. The solid line is a parametrization of Cronin effect derived from 
d+Au data in [22], the dotted lines reflecting the associated errors. 

momentum. 



3 J/ip mean transverse momentum 



Recombination is expected to modify transverse momentum distributions. To 
properly predict the modified spectra, one first needs to quantify the 
broadening coming from normal nuclear effects (Cronin effect). This effect was 
clearly seen at SPS by comparing p+p and p+A (p?), as well as in PHENIX at 
forward 3 rapidity [6]. At SPS, a simple parametrization could reproduce the 
(p?) values from p+p up to Pb+Pb: {pt)aa = (Pt)pp~^P°'^((Pt)) x ^ j where L is 
the average thickness of nuclear matter seen by a J/ip. The factor pa<5((p\)) 
stands for the nuclear density p, times the elastic gluon-nucieon scattering 
cross section a, times the average px kick given at each scattering 5((p^)). 
A review of SPS (including some FNAL data) was made in [22] together 
with a similar fit to RHIC forward data. This is presented as a solid line 
(with associated errors) on figure 4 (right). The shaded bands are predictions 
from [21], corresponding to either J/^'s from recombination (lower band) or 
to directly produced J/^'s (upper band). No clear sign for modification is seen 
and the A+A S((p"^)) can be interpreted in terms of normal broadening. 



3 The midrapidity p+p has too poor statistics to claim for a modification with 
respect to the d+Au measurement. 
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The rapidity spectra are also expected to be modified by recombination, but 
there is also no sign of this so far [10]. 

As a conclusion, I stress that RHIC and SPS data are not so easy to compare, 
even if they exhibit similar suppression at their highest energy densities. The 
amount of normal nuclear suppression is poorly known, especially at RHIC 
where it demands more d+A data. Nevertheless, RHIC preliminary suppres- 
sion seems anomalous and all the different models that can accommodate it 
suppose the formation of a QGP. To distinguish between them, a better pre- 
cision on data, and in particular on the kinematical distributions, is required. 
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